Gamma Rays Prom The Galactic Center and the WMAP Haze 
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Recently, an analysis ol data from the Fermi Gamma Ray Space Telescope has revealed a flux of 
gamma rays concentrated around the inner ~0.5° of the Milky Way, with a spectrum that is sharply 
peaked at 2-4 GeV. If interpreted as the products of annihilating dark matter, this signal implies that 
the dark matter consists of particles with a mass between 7.3 and 9.2 GeV annihilating primarily 
to charged leptons. This mass range is very similar to that required to accommodate the signals 
reported by CoGeNT and DAMA/LIBRA. In addition to gamma rays, the dark matter is predicted 
to produce energetic electrons and positrons in the Inner Galaxy, which emit synchrotron photons 
as a result of their interaction with the galactic magnetic field. In this letter, we calculate the flux 
and spectrum of this synchrotron emission assuming that the gamma rays from the Galactic Center 
originate from dark matter, and compare the results to measurements from the WMAP satellite. 
We find that a sizable flux of hard synchrotron emission is predicted in this scenario, and that this 
can easily account for the observed intensity, spectrum, and morphology of the "WMAP Haze" . 
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Dark matter in the form of weakly interacting massive 
particles (WIMPs) is predicted in many models to pro- 
duce observable signals such as gamma-rays, neutrinos, 
charged particles and synchrotron radiation through its 
annihilations pQ. Recently, an analysis of the first two 
years of data from the Fermi Gamma Ray Space Tele- 
scope has identified a component of gamma-ray emission 
that is highly concentrated within 0.5° around the galac- 
tic center, and which is spectrally and morphologically 
distinct from all known backgrounds 2J. While a new 
astrophysical source such as an undetected population of 
millisecond pulsars may create such a signal [3] , the mor- 
phology of this emission can most easily be accounted for 
by dark matter with a mass of 7.3-9.2 GeV annihilating 
primarily to tau leptons (possibly in addition to other 
leptons), and distributed in a cusped, and possibly adia- 
batically contracted [see e.g. |4], profile (with p oc r~ 1,18 
topcx r- 133 ). The normalization of this observed signal 
requires the dark matter annihilation cross section (to 
t+t~) to fall within the range of (av) = 4.3 x 10~ 27 to 
3.9 x 10 -26 cm 3 /s, similar to that required of a thermal 
relic [2J [5] . The mass range indicated by this observation 
is similar to that required to accommodate the excess 
events reported by the CoGeNT jS] collaboration, and the 
annual modulation observed by DAMA/LIBRA [8]. 
Recently, Buckley et al. [5] have shown that the phe- 
nomenology of this leptophilic dark matter model is con- 
sistent with known constraints from the relic abundance, 
collider physics and direct dark matter detection. 

If leptophilic dark matter is responsible for the gamma 
ray signal observed in the galactic center, then a sizable 
flux of energetic electrons/positrons will necessarily be 
injected into the inner galaxy. Such particles lose energy 
via synchrotron, inverse Compton and bremsstrahlung. 
Given that the dark matter's mass, profile, annihilation 
cross section, and annihilation channels are highly con- 



strained by the intensity, morphology, and spectral shape 
of the observed gamma rays from the galactic center, it 
is possible to make a robust prediction for the character- 
istics of the corresponding synchrotron signal, providing 
us the opportunity to either falsify or further support a 
dark matter interpretation of the gamma ray signal. 

GeV-scale cosmic ray electrons/positrons propagating 
in 10-100^G magnetic fields produce synchrotron ra- 
diation that peaks at GHz frequencies. The emission 
is thus naturally expected to fall within the frequency 
range studied by cosmic microwave background (CMB) 
missions, such as the Wilkinson Microwave Anisotropy 
Probe (WMAP) [9]. Data from WMAP and other CMB 
experiments can, therefore, be used to potentially con- 
strain or detect the occurrence of dark matter annihila- 
tions in our galaxy. Interestingly, WMAP observations 
have revealed an excess of microwave emission in the in- 
ner 20° around the center of the Milky Way, distributed 
with approximate radial symmetry [10] . Although this 
excess is controversial [11], and may have an intensity 
which is affected by the details of astrophysical fore- 
ground subtraction [12] . the spectrum and morphology 
of this source does not correlate with any known fore- 
ground |13J . Possible astrophysical origins, such as ther- 
mal bremsstrahlung (free-free emission) from hot gas, 
thermal dust, spinning dust, and Galactic synchrotron 
traced by low-frequency surveys have proven problem- 
atic [lOl Q2]. Alternatively, it has been suggested that 
the haze could be synchrotron emission from a distinct 
population of electrons and positrons produced by dark 
matter particles annihilating in the inner galaxy [141 115) . 

The possibility that light (^5-10 GeV) annihilating 
dark matter particles may explain the presence of the 
WMAP haze has been bolstered by recent claims that 
the magnetic field near the Galactic Center may be sig- 
nificantly stronger than previously thought. In partic- 
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ular, Crocker et al. [TB] analyzed the spectral break in 
the Inner Galaxy's non-thermal radio spectrum and set 
a lower limit of 50/iG on the magnetic field within 400 pc 
of the Galactic Center. Further work [T7| found a best fit 
for 100-300/iG magnetic fields in the region around the 
Galactic Center. Such large magnetic fields in the In- 
ner Galaxy would result in a brighter flux of synchrotron 
emission and a harder synchrotron spectral index than 
had previously been expected. 

In this study, we will address the question of whether 
the observed synchrotron emission from the inner kilo- 
parsecs of the Milky Way is consistent with the a dark 
matter interpretation of the gamma ray signal from the 
Galactic Center. The galactic dark matter distribution 
is thought to be consistent with an NFW profile [18], 
and in the case of the WMAP haze, it has been shown 
that small deviations from this distribution (such as an 
Einasto or Via Lactea profile), have little effect on the 
morphology of the WMAP haze [19]. However, in order 
to match the morphology of the 7-ray signal near the 
galactic center, we require a dark matter density profile 
which is steeper than a traditional r^ 1 employed in |18j . 
We note that such a contraction (to as steep as r -15 ) 
is predicted by models of adiabatic contraction due to 
baryonic interactions in the inner galaxy [3]. We note, 
however, that the position of the spectral break between 
the adiabatic contracted inner galaxy and the r _1 falloff 
of nearer the solar position is uncertain, and thus we con- 
sider two different dark matter halo profiles to bracket 
this uncertainty. First, we adopt a simple extrapolation 
of the distribution implied by the 7-rays from the galactic 
center (halo profile A): 

PA (r)=D GcY/cm^ 8 ^y, (1) 

where r is the distance from the Galactic Center, and we 
evaluate a = {1.18, 1.34}. The normalization value Do 
has been set such that the total mass within the solar 
circle is equal to that obtained with a r -1 profile and a 
local density of 0.3 GeV/cm 3 . For the value a = 1.18 we 
have Do = 0.273 GeV/cm 3 , and for a = 1.33 we obtain 
Do = 0.251 GeV/cm 3 . Additionally, we consider a broken 
power-law profile (halo profile B): 

PB (r) = 0.30 GeV/cm 3 ^ 8 " 5 ^ PC ^ , r> 0.5 kpc 

PB {r) = 5.1 GeV/cm 3 (^^j , r < 0.5 kpc. (2) 

We do not comment here on the case of the a = 1.18 
profile with a broken power law profile, as the results 
are not greatly altered compared to the unbroken pro- 
file. To generate the observed gamma ray flux from the 
Galactic Center, profile A requires an annihilation cross- 
section of 4.3 x 10~ 27 cm 3 s _1 for the case a=1.33 and 
3.2 x 10~ 26 cm 3 s _1 for the case a=1.18. For profile B, 



we evaluate only the case a=1.33 and set an annihilation 
cross-section of 1.3 x 10 -26 cm 3 s _1 [2]. 

It is possible that in addition to t + t~, the dark mat- 
ter may also annihilate to final states which produce few 
gamma rays, but contribute significantly to the cosmic 
ray electron/positron spectrum, and therefore to the re- 
sulting synchrotron emission. We will consider two cases: 
one in which the dark matter annihilates uniquely to 
t + t~, and another in which it annihilates equally to 
e + e~, and t + t~ (which we will refer to as "demo- 

cratic leptons"). The latter case injects roughly an order 
of magnitude more energy into cosmic ray electrons, and 
with a somewhat harder spectrum than that from taus 
alone. 

To model the propagation of electron and positrons 
produced through dark matter annihilation, we employ 
the cosmic ray propagation code Galprop |20l 121) , which 
calculates the synchrotron energy spectrum including ef- 
fects such as diffusion, reacceleration, and alternative en- 
ergy loss mechanisms including inverse Compton scatter- 
ing and bremsstrahlung radiation. We have modified the 
Galprop program to accept a dark matter annihilation 
spectrum from the DarkSUSY package [55], which has 
been formulated to produce the lepton fluxes for a variety 
of WIMP annihilation channels. Unless otherwise noted, 
we have adopted parameters identical to the best fit prop- 
agation model (base) determined in Ref . []15] , which were 
found to provide the best agreement to observed cosmic 
ray primary to secondary ratios. We adopt a diffusion 
constant D = 5.0 x 10 28 cm 2 s _1 , an Alfven velocity 
v Q = 25 km s _1 , and assume negligible convection. 

We parametrize the Galactic Magnetic Field according 
to: 

B(r,z) =Boe -(^-ikpc)/fl Se -(| z |-2kpc)/ z% ( 3 j 

where r and z represent the distance from the Galac- 
tic Center along and perpendicular to the the Galactic 
Plane. This exponential distribution follows the standard 
set forth in [50] [5T] , and we employ values of R s and z s 
in order to approximately match the magnetic field dis- 
tribution both near the galactic center set by |16j . and 
locally set by [53]. We note that these parameters suffer 
from significant uncertainties, and we demonstrate the 
impact of these uncertainties on our models below. 

In order to extract the WMAP haze residual from the 
WMAP dataset, we utilize the foreground subtraction 
templates employed in Ref. [13 using the template sub- 
traction method CMB5, and we complete the mean sub- 
traction of the background dataset as in Ref. [19] . Since 
the dominant source of error in the determination of the 
WMAP haze stems from the accuracy of these subtrac- 
tion methods - rather than from the measured error in 
the WMAP dataset, we resist calculating error bars for 
our measurement, as they would serve primarily to belie 
to actual error in the WMAP haze measurement (for a 
quantifications of these errors, see Fig. 8 in [13]). In- 
stead, we include one- a error bars based on the temper- 
ature fluctuations in each radial and energy bin. In order 
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FIG. 1: Synchrotron emission from dark matter annihila- 
tions as a function of latitude below the Galactic Center for 
8 GeV dark matter particles annihilating to t + t~ , and us- 
ing halo profile A with a slope a = 1.18. The total cross 
section of av = 3.2 x 10 -26 cm 3 /s was chosen to normal- 
ize the gamma ray signal observed from the Galactic Center. 
The magnetic field model used is given by B(r,z) = 27.7 jj,G 
e _ (r _ lkpc)/3 . 2kpc e _(|»|_2kpo)/i.akpc < The " dark matter " sig _ 

nal corresponds to the synchrotron signal provided by dark 
matter leptons in our simulations, while the total signal in- 
cludes the "dark matter" signal added to the zero-point offset. 
The error bars shown correspond to the 1-cr temperature fluc- 
tuations in each radial and energy bin [13] 



to model the isotropic component which was subtracted 
from the WMAP dataset in the template fits, we add an 
isotropic component to the simulated haze which is al- 
lowed to float individually in each energy band (labeled 
the zero point offset). This is intended to account for 
systematic uncertainties in the template subtractions, as 
well as for any other astrophysical sources of hard syn- 
chrotron that might be present. 

Our models indicate that the details of the annihilation 
pathway postulated by [2] are the most important fac- 
tor in controlling the dark matter match to the WMAP 
haze. We find that the r-only annihilation channel is 
not able to reproduce the intensity of the WMAP haze 
using a = 1.33 with either profile- even if we allow for 
an unrealistically large magnetic field of 180^tG in the 
haze region. However, In the case of an a = 1.18, we are 
assisted by the much larger cross-section of 3.2 x 10~ 26 
cm s _1 required to match the observed 7-ray signal. In 
Figure [l] we show that this allows for a reasonable match 
to the WMAP haze if we allow for strong magnetic fields 
on the order of 28/iG in the haze region. We find that 
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FIG. 2: As in Fig. [l] for a dark matter model which an- 
nihilates democratically into leptons with a cross section of 
(7v = 9.6 x 10 -26 cm 3 s _1 (again, normalized to the observed 
gamma ray signal). The magnetic field model used is given 
by B(r,z) = 6.28 /iG e -(—ikpc)/4.0k P c e -(|*|-2k P c)/2.2kpc^ 



the results are changed only slightly if we introduce a 
spectral break for this profile. 

In the case where dark matter instead annihilates 
democratically into leptons, we obtain a strong enhance- 
ment to the intensity of the WMAP haze for two reasons: 
(1) the cross-sections are multiplied by three because we 
obtain negligible direct 7-ray production from the muon 
and electron channels, so it is only the r cross section 
which must be held constant in [2], (2) the electron chan- 
nel converts energy entirely into leptons as opposed to the 
much smaller contributions from the r and fi channels. 
In order to account for this much larger input lepton 
density, we must greatly decrease the magnetic field to 
maintain a constant synchrotron intensity. In Figure [2] 
we show the match for profile A with a — 1.18, which 
employs a magnetic field of strength 6.3/iG in the Haze 
region. 

For the case of a = 1.33, we can observe the small 
differences between profile A and profile B. In Figure [3] 
we show the best fit dark matter profile for profile A with 
a cross-section of 1.3 x 10 -26 cm 3 s _1 , using a magnetic 
field of strength 9.4/xG in the haze region. This contrasts 
only slightly with our result in Figure |4j where we employ 
profile B and a cross-section of 3.9 x 10 -26 cm 3 s _1 . Thus 
we find the effect of both extremes of our dark matter 
profile extrapolation to have only a minimal effect. 

In summary, a spectrum of gamma rays sharply peaked 
at 2-4 GeV from the inner 0.5° around the Galactic Cen- 
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FIG. 3: As in Fig. [2] for a dark matter profile with 
a — 1.33 and profile A. The total annihilation cross- 
section has been normalized to av — 1.3 x 10 -26 cm 3 s -1 
The magnetic field model is given by B(r,z) = 9.43 fiG 
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FIG. 4: As in Figs. [3] except with profile B. The total annihi- 
lation cross-section has been normalized to av = 3.9 x 1CP 26 
cm 3 s _1 The magnetic field model is given by B(r,z) = 9.41 /xG 

e _( r ._lkpc)/4.0kpc e -(|z|-2kpc)/2.2 

ter has recently been identified within the data of the 



Fermi Gamma Ray Space Telescope [2J. If these gamma 
rays are interpreted as dark matter annihilation prod- 
ucts, this tightly constrains the dark matter particle's 
mass, annihilation cross section, annihilation channels, 
and halo profile. In this letter, we calculate the syn- 
chrotron emission from the Inner Galaxy predicted in this 
scenario, and compare it to the observed features of the 
WMAP Haze. Although the fits we have shown are not 
perfect matches to the data, they are well within the ex- 
pectation given the likely imperfections in simulations of 
the dark matter profiles, and lepton diffusion employed. 
We further note that the systematic uncertainties in the 
determination of the WMAP haze signal exceeds the er- 
ror bars plotted here, and possibly introduces morpholog- 
ical and spectral uncertainties [13] . We find agreement 
with the observed intensity, spectrum, and morphology 
of WMAP haze, especially if dark matter is seen to an- 
nihilate democratically to leptonic final states (equally 
to e + e~, and t + t~). If dark matter annihila- 

tions produce only t + t~ , we can still match the observed 
WMAP intensity if we resort to relatively strong mag- 
netic fields in the haze region. We note the magnetic 
field strengths between these two scenarios differ greatly, 
and an accurate determination of the magnetic field in 
this region would constrain the possible dark matter an- 
nihilation pathway. 

We emphasize that these fits to the characteristics of 
the WMAP Haze were obtained with relatively little free- 
dom in the astrophysical or dark matter parameters. In 
particular, the mass, annihilation cross section, and halo 
profile are each tightly constrained by the observed fea- 
tures of the Galactic Center gamma ray signal. We note 
that we have allowed the magnetic field energy density 
to modulate in order to obtain best fits to the WMAP 
dataset. Although these magnetic field choices allowed 
us to adjust the morphology and spectrum of the of the 
synchrotron emission to a limited degree, we had only 
moderate ability to adjust the overall synchrotron inten- 
sity. Finally, we note that these magnetic field modu- 
lations are observationally testable. For this reason, we 
find it particularly interesting that the dark matter model 
implied by the observed gamma ray signal so naturally 
yields synchrotron emission consistent with the WMAP 
Haze. While these results our not unique to our dark 
matter models as compared to those of more massive par- 
ticles (see e.g. [HI El]), these results show the promise 
for radio observations to constrain the scenarios of high 
energy dark matter annihilation such as those shown in 
[2], and the need for further studies, such as a model for 
synchrotron emission near the galactic center where the 
total 7-ray intensity is much better constrained. 
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